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Abstract 

We discuss short-baseline electron and muon neutrino disappearance searches into sterile 
neutrinos at a Very Low Energy Neutrino Factory (VLENF) with a muon energy between 
about two and four GeV. A lesson learned from reactor experiments, such as Double Chooz 
and Daya Bay, is to use near and far detectors with identical technologies to reduce the sys- 
tematical errors. We therefore derive the physics results from a combined near-far detector 
fit and illustrate that uncertainties on cross sections x efficiencies can be eliminated in a 
self-consistent way. We also include the geometry of the setup, i.e., the extension of the 
decay straight and the muon decay kinematics relevant at the near detector, and we demon- 
strate that these affect the sensitivities for Am 2 > 30 eV 2 , where oscillations take place 
already in the near detector. Compared to appearance searches, we find that the sensitivity 
depends on the locations of both detectors and the muon energy, where the near detector 
should be as close as possible to the source, and the far detector at about 500 to 800 m. In 
order to exclude the currently preferred parameter region, at least 10 19 useful muon decays 
per polarity are needed for = 2 GeV, or, alternatively, a higher muon energy can be 
used. 
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1 Introduction 



The test of three-flavor neutrino oscillations in solar, atmospheric, long-baseline, and reactor 
experiments has been so far very successful, where a non-zero mixing angle #13 has been re- 
cently established by the Daya Bay and RENO reactor experiments above the 5a confidence 
level [TJ[2] . On the other hand, neutrino oscillations at short baselines, i.e., L <C E/Am^ 
where atmospheric oscillations have not yet developed, face a tension between the strong 
constraints from many short-baseline experiments, and several observed anomalies which 
may be described by eV-scale sterile neutrinos. In particular, short-baseline electron neu- 
trino disappearance may cause an anomaly identified in Gallium experiments (3|, electron 
antineutrino neutrino disappearance may lead to lower than predicted reactor antineutrino 
fluxes pip)] , and electron antineutrino appearance may be driven by sterile neutrinos in the 
LSND |6j and MiniBooNE [7] experiments. In the simplest models, one would add one 
extra sterile generation to fit these data by a 3 + 1 model, separated by Am 2 ~ leV 2 , 
see, e.g., Ref. [8). However, it turns out that the tension between MiniBooNE neutrino 
and antineutrino data (CP violation cannot be described by this model), and the tension 
between appearance and disappearance data basically rule out this model, see, e.g., Ref. [9]. 
Therefore, 3 + 2 models have been more recently used |10|, which allow for CP violation 
and include additional degrees of freedom. In any of the above models, crucial information 
comes from the disappearance channels, which may be the "cleanest" channels to measure 
the oscillation parameters (see Sec.[2]for details). In addition, electron neutrino and antineu- 
trino disappearance searches are needed to directly test the Gallium and reactor anomalies, 
respectively. 

Various new experiments have been proposed to test short-baseline neutrino oscillations, 



see Appendix A of Ref. 11 for a recent summary of alternatives. In this study, we focus 



on a very low energy neutrino factory (VLENF), which is a neutrino factory with a low 
muon energy of about two to four GeV which does not require muon cooling or muon accel- 



eration 12 , and could be the first phase of a staged neutrino factory program. Compared 



to muon neutrino appearance, discussed in Ref. 12 , we discuss the electron and muon 
neutrino (antineutrino) disappearance channels. These are qualitatively different from the 
appearance channels because of different systematics: while the appearance channels are 
limited by (charge mis-identification and neutral current) backgrounds, the disappearance 
channels suffer from unknown cross sections x efficiencies. Therefore, similar to reactor 
experiments such as Double Chooz or Daya Bay, near and far detectors have been proposed 
for the high energy neutrino factory to control the systematical errors [13] . In addition, 



near detectors very close to the muon storage ring experience geometry effects from the 



extension of the decay straight and beam divergence 14 . In Ref. |13] the geometry effects 
from the decay straight have been taken into account, which lead to a smearing of the 
oscillation probabilities, whereas the detectors were assumed to be far enough away from 
the source (or small enough) not to experience any detector geometry effects (far distance 
limit). However, for the substantially lower muon energy of the VLENF, the beam will be 
larger from the muon decay kinematics only, which means that the detector geometry has to 
be taken into account. In addition, the near and far detectors are supposed to be as similar 
as possible, which is difficult to reconcile with different detector diameters. Therefore, we 
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will simultaneously integrate over straight geometry and detector surface area in this study, 
and show the impact of these effects. We discuss electron neutrino disappearance for most 
of this study, since it is directly relevant to test some of the anomalies, and we point out 
the differences for muon neutrino disappearance at the end. 

This study is organized as follows: we recapitulate the phenomenology of sterile neutrino 
disappearance searches in Sec. |2j Then in Sec. [3j we introduce our setup, systematics 
treatment, and geometry treatment. In Sec. |4j we illustrate the impact of the beam and 
detector geometry, and of the systematics assumptions. Then in Sec. [5j we perform a two- 
baseline optimization of near and far detectors. We show our main results for electron and 
muon neutrino disappearance in Sec. [6j and we conclude in Sec. [7} 

2 Sterile neutrino phenomenology 

In most experiments, the disappearance or appearance of neutrinos at short distances is 
described in the two-flavor limit 

P aa = l-sin 2 (2# cff )sin 2 A, (1) 
P a p = sin 2 (2# eff )sin 2 A, (2) 

where A = Am 2 L/(4:E) and 6 e s is an effective mixing angle. However, this description 
cannot be used for a self-consistent description of a multi-channel experiment. In order to 
demonstrate that, consider the simplest possible case, a 3 + 1 framework with one extra 
sterile neutrino. Although this case is basically excluded, it is useful two illustrate some 
of the considerations to be taken into account here which also apply to more sophisticated 
models. The parameterization-independent probabilities in the limit A41 3> A31 ~ (A^- = 



AmjjL/ (4E)) can be written as 

P ee = 1 - 4|t/ e4 | 2 (l - |C/ e4 | 2 ) sin 2 A 4 i , (3) 

P w = l-4|fV| 2 (l-|fV| 2 )sin 2 A 41 , (4) 

Pe» = P^e = 4|[/ e4 | 2 |[/ At4 | 2 sin 2 A 41 , (5) 



where we show the most interesting channels for the Neutrino Factory. From this simple 
model, it is immediately clear that LSND-motivated electron or muon neutrino appearance, 
which requires |C/ e4 | > and \U^\ > 0, must be accompanied by electron and muon neutrino 
disappearance, and that the disappearance searches provide important and strict constraints 
on theoretical models. Especially, the disappearance probabilities are proportional to jf/^j 2 , 
the appearance probabilities to |t/o, 4 | 4 , i.e., the appearance probabilities are suppressed by 
two more powers of the new mixing angles. In a specific parameterization, electron and 
muon disappearance can be described by the two-flavor limit Eq. ([!]) by different mixing 
angles (such as 9u and 6* 24 , respectively), whereas the appearance probability in Eq. (|2| is 
proportional to the product of these mixing angles; see e.g. Ref. [15] for a direct comparison. 
This means that Eq. Q or Eq. @ can be used independently to effectively describe an 
individual disappearance or appearance channel. If, however, the information from different 
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channels is to be combined, there will be a model-dependent (but well defined) interplay 
among the channels. 

Now what are the consequences for the Neutrino Factory? For electron or muon neutrino 
appearance, the main limitations are charge mis-identificatiorj^] and neutral current back- 
grounds. At least for small mixing angles, other systematics, such as the flux and cross 
section uncertainties, are less relevant. That is quite fortunate, since it would be probably 
hard to quantify this systematics. Consider, for instance, a near detector to measure the 
cross sections for muon neutrinos. For large enough Am| 1; oscillations may already take 
place in the near detector. Unless it is clearly defined how the theoretical interplay between 
muon neutrino disappearance at the near detector and muon neutrino appearance at the 
far detector works, one cannot disentangle the cross sections from oscillation physics in that 
case. Thus, an effective two-flavor description is not sufficient to obtain reliable information 
on this systematics. 

For the disappearance channels, the situation is very different. Here the electron and muon 
neutrino disappearance can be described by Eq. ([I]) in the effective two-flavor limit with 
different effective mixing angles. Because small deviations from unity are to be measured, 
a lesson learned from reactor experiments 16 , 17] , such as Double Chooz or Daya Bay, has 



been to use (more or less) identical near and far detectors to cancel systematical errors. 
Compared to reactor experiments, where the flux is the unknown, the cross section x ef- 
ficiency uncertainties are the dominant systematics to be canceled by the near detector. 
From the oscillation framework point of view, since both near and far detectors measure the 
same flavors, the oscillation probabilities in both detectors can be described by the same 
probabilities. In the 3 + 1 model, one can simply use Eq. M for v e or disappearance. 
In a 3 + iV model, the translation into the model parameters is more complicated. Nev- 
ertheless, it is clear that any detected difference between near and far rates (or neutrino 
and antineutrino rates) must be due to oscillations, and cannot come from unknown cross 
sections or efficiencies. Very interestingly, the near and far detectors may change their roles 
as a function of Am^. While for small Am^ the far detector measures the oscillations and 
the near detector the normalization, for large Am^ the near detector measures the oscilla- 
tions and the far detector, where the oscillation averages out, the normalization [13). We 
choose an effective two-flavor framework in the following to quantify the performance, and 
we assume CPT invariance for the sake of simplicity (the disappearance channels are always 
CP invariant). Note that nevertheless CPT invariance tests of the disappearance channels 



are well motivated, see Ref. 13 for the Neutrino Factory. We do not perform a combined fit 
of appearance and disappearance channels, or electron and muon disappearance channels, 
since such a fit can only be performed for a specific model. Note, however, that the final 
experiment can of course be used to test more complicated scenarios. 
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Figure 1: Possible geometry of the VLENF near-far setup for = 2 GeV (not to scale). 
3 Setup and simulation techniques 

For the simulation, we choose a near-far detector setup with a near detector of 200 1 fiducial 
mass x efficiency, and a far detector of 1 kt fiducial mass x efficiency, see Fig. [T] for a 
possible geometry. Since the detector geometry will be important for close distances, we 
assume cylindrical shapes with diameters of 6 m (perpendicular to the beam axis). We test 
two different muon energies: = 2 GeV and E^ = 4 GeV. The decay straight length s is 
assumed to be 100 m for E^ = 2 GeV, and 200 m for E^ = 4 GeV. For very short baselines 



and line neutrino sources, the baseline L is ill-defined 14 . Therefore, the detector locations 
are specified by the distance d to the end of the decay straight, i.e., d + s < L < d. For 
the integrated luminosity, we use 10 19 useful muon decays per polarity, and 10 18 useful 
muon decays in some cases where explicitely specified. For the detection threshold, we 
use 0.25 GeV and 0.5 GeV for = 2 GeV and E^ = 4 GeV, respectively. For the energy 
resolution of the detectors, we use AE/GeV = 0.1 y/E/GeV for u e (u e ) disappearance, and 
AE/GeV = 0.05 ^/E/GeV for Up (v^ disappearance, which is motivated by a totally active 
scintillator or a liquid argon detector. Neutral current backgrounds are included at the level 
of 10 -4 , which have, however, only a very small effect. Note that a magnetization of the 
detector for the disappearance searches may or may not be necessary, depending on the 
underlying physics model, i.e., what happens to the other flavor. For instance, in the 3 + 1 
model discussed above, one may have large z/ M disappearance driven by non-zero \U^\, see 
Eq. Q, whereas appearance vanishes at the same time for |{7 e4 | = 0, see Eq. Since 
we only use the disappearance channels, we assume that charge mis-identification is either 
under control by a magnetic field, or suppressed by the underlying physics. Note that the 
chosen parameters are consistent with the currently discussed ones of the VLENF study 



group [18 ; see also Ref. 12 



For the geometric treatment of the neutrino line source and detector geometry, we follow 
Refs. |13|[l4], and for the simulation, we use the GLoBES (General Long Baseline Experi- 



ment Simulator) software 19 20 . GLoBES (up to version 3.1) assumes that the detectors 



^or instance, for neutrino production by /i + — > e + + v e + v^, one has to distinguish a few v e — > 
charged current events from the dominant — > event rate by charge-identification of the produced 
muons (anti-muons) by means of a magnetic field. 
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are far enough away from the source to treat the source as point source, and that the detec- 
tors are small compared to the beam divergence. If we start from the differential event rate 
from a point source dNps/dE without oscillations, as the one used in GLoBES, we can take 
into account the extension of the straight and the detector and the effect of oscillations by 
an averaged event rate 

d+s d+s 

dN avg _l f dN 1 f dN PS (L,E) 



dE s 



f dN 1 f dN PS (L,E) 

/ dE s / dE £a(L,E)P aa (L } E)dL. (6) 



Here e a (L,E) = A e fi ia /An et parameterizes the integration over the detector geometry for 
a fixed baseline L, energy E, and flavor z/ Q , where A eS ^ a (L, E) is the (energy and flavor 
dependent) effective area of the detector, and A Det is the actual surface area. Note that 
the oscillation probability P aa appears inside the integral, since different parts of the decay 
straight contribute differently. In addition, note that it is assumed that the differential muon 
decay rate per (straight) length is equal along the decay straight. Since dN P $/dE oc l/L 2 , 
we can re-write this as 

dN avg _ dN PS {L eS ,E) {L, E) _ dN PS (L eS ,E) 



dE dE s 



J^P m (L,E )d L= iN *^ P ( E) (7) 



d 

with the average efficiency ratio times probability 

d+s 



P{E) = *f f P aa {L, E)dL (8) 



and the effective baseline 

L eS = ^d(d + s) , (9) 

which is defined such that P(E) = 1 for e a (L,E) = 1 and P aa (L,E) = 1. As a conse- 
quence, one would use the usual detector definitions in GLoBES with the effective baseline 
L e ff, which are to be corrected by Eq. We compute Eq. ^ directly in a user-defined 
probability engine in GLoBES, including both neutrino source and detector geometry. For 
the sake of simplicity, we assume that the (machine-dependent) beam divergence is smaller 
the the beam spread given by the muon decay kinematics. In this ideal case, one can easily 



compute e a (L,E) independent of machine-dependent effects more or less analytically 14 
We will demonstrate that there is a substantial effect coming from the extension of the 
beam compared to the detector, whereas machine-dependent effects may lead to additional 
smearing if the divergence is not under control. Note that the intrinsic effect of the muon 
decay kinematics cannot be removed, even in an ideal machine. 

For the systematics treatment, we follow the reactor experiments with two (or more) detec- 



tor, rather than the usual neutrino factory description; see Refs. 13 21 for details. From 



the reactor experiments, we known that short baseline electron neutrino disappearance is 



mostly affected by the signal normalization uncertainty (see, e.g., Refs. 17 21 for reactor 



experiments). Here, compared to the reactor experiments, our signal normalization error is 
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not dominated by the flux, which may be known at the level of 0.1% using various mean 
monitoring devices [22], but the knowledge of the cross sections x efficiencies. Because our 
neutrino energies span the cross section regimes from quasi-elastic scattering, over resonant 
pion production, to deep inelastic scattering, it is difficult to estimate the degree the cross 
sections will be known at the time of the measurement. The efficiencies depend on the 
detection processes and detector properties, which means that their uncertainties are also 
difficult to pin down. For reactor experiments, on the other hand, the inverse beta decay 
cross sections are well known, but the fluxes are relatively uncertain. Both classes of ex- 
periments have, however, in common that these uncertainties can be controlled by using 
detectors as identical as possible. In our case, where the near and far detector masses are 
different, the far detector may consist of five modules identical to the near detector, as it 
is illustrated in Fig. [I] We neglect the extension of the detector along the beam axis, since 
it is expected to me much smaller than the extension of the decay straight. However, in 
practice, the location of the interaction vertex can be measured to some degree. 

We adopt the most conservative case for systematics, which is that the cross sections x 
efficiencies (x flux) are fully uncorrelated among the bins, but fully correlated between the 
near and far detectors]^] This assumption is conservative because is corresponds to cross 
sections x efficiencies with an unknown shape error, where the shape is to be measured by 
the near detector. Unless noted otherwise, we assume that cross sections x efficiencies are 
only known to the level of 10% (within each bin), where even larger errors do not matter 
in the oscillation region. In addition, we use a normalization error uncorrelated between 
near and far detectors, but fully correlated among the bins, which may come from a fiducial 
mass uncertainty. This error is, for reactor experiments, believed to be controlled below the 
percent level. We use 0.6% [21] , whereas Daya Bay claim that they can even do significantly 
better (0.2% [I]). Since it depends on the detector properties and detection interactions, 
a more conservative choice seems reasonable. Finally, the backgrounds are assumed to be 
known within 35% and the energy calibration to 0.5%, uncorrelated between the detectors. 
In the next section, we will discuss where these uncertainties are important. 

4 Impact of geometry and systematics 

Let us now study the impact of geometry and systematics, where we use = 2 GeV, 
d\ = 20m (200 1), and d 2 = 500m (lkt). For the most of the following, we focus on v e 
ip e ) disappearance, and discuss the differences to muon neutrino disappearance at the end. 
In Fig. [2J left panel, the exclusion region in sin 2 29 and Am 2 (right hand sides of curves) 
is shown for different geometry assumptions. As we will demonstrate below, this setup 
is optimized for not too small values of Am 2 . The curve "no systematics" represents the 
statistics limit of the far detector only, without systematics or backgrounds. It is simulated 
in the point source and far distance approximations, which means that the baseline is 
computed with Eq. ^ and that e a (L,E) = 1 in Eq. ([8]), respectively, as it is usually done 
for far detectors. The curve "point source" is obtained for the same geometry assumptions, 

2 We use 17 bins with bin widths following the energy resolution of the detector. To avoid aliasing effects, 
we use the built-in filter of GLoBES on the oscillation probability. 
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Figure 2: Exclusion region in sin 2 29- Am 2 (right hand sides of curves) for v e disappearance for different 
geometry assumptions (left panel) and systematics assumptions (right panel); see main text for details 
(90% CL, 2 d.o.f.). The curves "no systematics" represents a single detector at d = 500m using statistics 
only, whereas the other curves correspond to near-far detector setups, where the red thick curves include 
(conservative) full systematics and geometry effects. Here = 2 GeV, 10 19 useful muon decays per polarity, 
dt = 20 m (200 1) and d 2 = 500 m (lkt). 

but for a near-far detector simulation including full systematics. The first peak of the 
sensitivity at about Am 2 ~ 3 eV 2 corresponds to the far detector, as it can be easily seen. 
In this case, the near detector measures the normalization, the far detector the oscillation. 
The near detector sensitivity peaks at about Am 2 ~ 13 eV 2 , where, however, oscillations 
are still present in the far detector (see "no systematics" curve). Therefore, the optimal 
sensitivity is reached at somewhat larger values of Am 2 , where oscillations in the far detector 
average out and the cross sections x efficiencies are safely measured in the far detector. In 
this case, near and far detectors swap their roles. This swapping is also the reason why it is 
difficult to obtain reliable sensitivity predictions for effective far detector simulations only, 
especially for Am 2 3> leV 2 , where oscillations take place in the near detector. As the next 
step in Fig. [2j left panel, we take into account the extension of the decay straight in the 
curve "straight averaged", which leads to some averaging in the region the near detector 
is most sensitive to. This case corresponds to using Eq. (Tsj) with e a (L,E) = 1. The final 
result, curve "straight and detector averaged", shows the additional impact of the detector 
geometry, i.e., the full Eq. In this case, the sensitivity for large Am 2 , coming from the 
near detector, almost vanishes. The reason is that the detected spectrum effectively peaks 
at lower energies due to the muon decay kinematics, which the near detector (but not the 
far detector) is sensitive to. For very large Am 2 3> 100 eV 2 , the oscillations average out in 
both near and far detectors, and the sensitivity is given by the externally imposed error of 
10%. 
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We have also tested if one can reduce the effect of the straight averaging by the use of two 
beam current transformers (BCTs), one before and one after the straight. One may assume 
that the number of muon decays along the straight, which is proportional to the difference 
between the two beam currents, is exponentially distributed along the straight. However, 
already at = 2 GeV, the mean lifetime (decay length) of the muon is about 6 km, which 
means that the effect is small and not visible in the sensitivities. 

In the right panel of Fig. |2j we show the impact of different systematics assumptions, where 
again the "no systematics" and full systematics ("10% shape error") curves are used for 
reference. First of all, we reduced the systematical errors one by one. As expected, a 
smaller shape error (curve "2% shape error") improves the sensitivity for large Am 2 , where 
the oscillations in near and far detectors are averaged out and this systematics dominates. 
None of the other systematics has an effect on the sensitivity if improved separately. Only 
if all of the systematical errors can be controlled better, the curve "low systematics" is 
obtained]^] Since it is difficult to say how realistic this case is, we rely on our standard (more 
conservative) values in the following, unless explicitely stated otherwise. In Fig. [2j we also 
show the result for the far detector only (gray dashed-dotted curve), from which one can 
see that the near detector is necessary to maintain the sensitivity for low Am 2 ~ eV 2 in the 
presence of the systematical errors. We include the decay straight and detector geometries 
in the following calculations, as well as the near detector, where the full systematics curve 
("10% shape error") represents our standard values. 

5 Two-baseline optimization 

Whereas the muon energy is limited by other constraints, such as accelerator and storage 
ring, the detector locations of the two detectors can be almost freely chosen. We therefore 
show in Fig. [3j left panel, the two-baseline optimization for v e disappearance. The different 
regions correspond to optimal detector locations for the depicted values of Am 2 , where the 
respective reaches in sin 2 29 are given in the figure caption. For Am 2 = leV 2 , two optimal 
regions are found, where two central choices are marked by points A and B. Point A is 
within a region where the near detector is as close as possible to the source, and the far 
detector in a distance between about 200 m and 1 km. Point B corresponds to a longer far 
detector baseline which helps for small values of Am 2 . However, a somewhat farther near 
detector distance is preferred, where the near detector also adds to the sensitivity directly. 
The larger Am 2 is, the smaller distances for the far detector are preferred, where in all cases 
the near detector may be as close as possible to the source. Point C is representative for 
a region with better sensitivity for Am 2 > 10 eV 2 . The asymmetry in Fig. [3] (left panel) 
with respect to the symmetry axis d± = d,2 comes from the different detector masses. For 
equal detector masses, we do not find a qualitative differences apart from the plot becoming 
symmetric with respect to this axis. 

The results in the Am 2 -sin 2 26 l -plane are shown for the marked setups in Fig. [3J right 
panel. Setup B has the best sensitivity for small values of Am 2 , but for larger Am 2 it is 

3 Low syst.: Fiducial volume/normalization error 0.1%, shape error 2%, calibration error 0.1%, back- 
ground error 10%. 
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Figure 3: Left panel: Optimal detector distances for v e disappearance, = 2GeV, 10 19 useful muon 
decays per polarity. The contours correspond to sensitivities in log 10 (sin 2 29) down to -1.9 (Am 2 = 1 eV 2 ), 
-1.95 (Am 2 = 10eV 2 ), -1.5 (Am 2 = 32eV 2 ), -1.2 (Am 2 = lOOeV 2 ) for the depicted values of Am 2 (90% 
CL, 2 d.o.f.). Right panel: Sensitivities in sin 2 29-Am 2 for the "optimal" setups marked in the left panel. 
For comparison, setup A is shown with 10 18 useful muon decays per polarity, and the best-fit region from 
Ref. [4] (Fig. 8) is shown as gray-shaded region (99% CL). 

not optimal. Setup C is best for Am 2 > leV 2 , as expected, whereas setup A is a good 
compromise between the small and large Am 2 sensitivities. Therefore, we have chosen it 
for reference. In order to check if the chosen setup match the needs for v e disappearance, 
we show the best-fit region from Ref. |2] (Fig. 8) as gray-shaded region for reference (99% 
CL). One can easily read off the figure that all setups can exclude this best-fit region very 
well, with setup C actually covering the largest part. Note again that the very large Am 2 
coverage is limited by the external knowledge on cross sections x efficiencies for the VLENF, 
whereas the flatness of the best-fit region for large Am 2 simply means that the reactor 
experiments cannot resolve the oscillations (necessary to exclude this part). Therefore, one 
should probably not overemphasize this part. For reference, we also show in Fig. [3j right 
panel, the curve for point A with 10 18 useful muon decays (per polarity) only, as different 
luminosities are currently being discussed. One can clearly read off the figure, that the 
statistics is not sufficient to fully exclude the interesting part of the best-fit region at lower 
values of Am 2 . 

In Fig. [4j we perform a similar analysis for Ep = 4 GeV. In the left panel (two-baseline 
optimization), only the the point A is chosen as for E^ = 2 GeV. The qualitative results 
of the two-baseline optimization are however similar to the above case, with the exception 
that the optimal region with a long baseline has disappearance in the chosen baseline win- 
dow. In addition, for very large Am 2 , somewhat longer baselines are preferred to avoid 
the geometry effects. In the right panel, the sensitivities for the chosen points are shown, 
with rather similar results. Again, point A appears to be a good compromise between the 
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Figure 4: Left panel: Optimal detector distances for v e disappearance, = 4 GeV, 10 19 useful muon 
decays per polarity. The contours correspond to sensitivities in log 10 (sin 2 29) down to -2.15 (Am 2 — 1 eV 2 ), 
-2.5 (Am 2 = 10 eV 2 ), -1.85 (Am 2 = 32 eV 2 ), -1.45 (Am 2 = 100 eV 2 ) for the depicted values of Am 2 (90% 
CL, 2 d.o.f.). Right panel: Sensitivities in sin 2 29-Am 2 for the "optimal" setups marked in the left panel. 
For comparison, setup D is shown with 10 18 useful muon decays per polarity, and the best-fit region from 
Ref. [4] (Fig. 8) is shown as gray-shaded region (99% CL). 



small and large Am 2 sensitivities, but point D performs better for small Am 2 . The ab- 
solute sensitivities are significantly better than in the = 2 GeV case, which is a result 



qualitatively different from the appearance optimization in Ref. 12 , because statistics are 
important for the disappearance channels. In this case, even the low luminosity curve with 
10 18 useful muon decays covers the relevant parameter space. In summary, for electron neu- 
trino disappearance, either 10 19 useful muon decays or = 4 GeV are required to exclude 
the discussed best-fit parameter space region. 



6 Results and comparison 



We summarize our results in Fig. [5] for disappearance (left panel) and v e disappearance 
(right panel). For disappearance (left panel), optimized setups for different values of 
are compared. As discussed above, the sensitivity for E^ = 4 GeV is significantly better than 
for Ep = 2 GeV, but both setups can in principle cover the relevant parameter region. The 
large Am 2 region coverage, marked "systematics limit", depends on the assumed external 
knowledge of cross sections x efficiencies. The modification of the E^ = 2 GeV case for 
an improved 2% error is shown as dashed curve in the same color. In this case, one can 
clearly fully cover the discussed (gray-shaded) parameter space. In addition, we show two 
curves for alternative approaches to u e disappearance measurements. One example is a low 



7 ~ 30 beta beam 24 , shown as dashed curve. This setup is in a way very similar to ours 
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Figure 5: Left panel: Exclusion region in sin 2 29- Am 2 (right hand sides of curves) for u e disappearance 
with 10 19 useful muon decays per polarity for = 2 GeV (point A, 20 m + 500 m) and — 4GeV 
(point D, 20 m + 800 m) at the 90% CL, 2 d.o.f. The orange-dashed curve shows the = 2 GeV result 
for an improved (2%) shape error. For comparison, two reference setups are shown: a radioactive ion 
facility [23] (dashed-dotted red curve in Fig. 6 therein; 10 14 8 Li ions per second) and a low gamma beta 
beam 24 (red curve in Figs. 4 and 5 therein). Right panel: Comparison between disappearance (solid) 
and v e disappearance (dotted) for the three test points (10 19 useful muon decays per polarity, — 2 GeV; 



90% CL, 2 d.o.f). The combined SciBooNE and MiniBooNE disappearance result from Ref. 25 is shown 
for comparison. 

both in terms of 7 (our 7 ~ 19 for = 2 GeV) and the beam geometry. For the detection 
reaction, however, inverse beta decay is used, and thus the systematical error is assumed to 
be controllable at the level of 1% (ours: 10%). Furthermore, this detection process is only 
sensitive to u e . Another example, shown as dashed curve, is the radioactive ion facility in 



Ref. 23 . Here 7 ~ 1 ions are injected into a An detector. Again, the systematical error is 
assumed to be 1%. The final result depends on the ion intensity, where we have shown the 
most aggressive scenario in Fig. [5] (10 14 ions per second). Here the main issue seems to be 
the coverage of the low Am 2 region, which requires relatively long baselines (which cannot 
be realized within such a detector). 

For Up disappearance (right panel, solid curves), we have used the same fiducial masses x 
efficiencies for the sake of simplicity. Therefore, the main difference to u e disappearance are 
the beam spectrum peaking at higher energies (for the same E^), and the better energy 
resolution. For the optimization, we have not find any qualitative differences compared to 
u e appearance, apart from the fact that slightly longer baselines are preferred, especially for 
point B in Fig. [3] For the sake of consistency, we therefore show the same optimization points 
in Fig. [5] as in Fig. |3l and we show the v e disappearance (dotted curves) for comparison. One 
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v e disappearance 



disappearance 




Figure 6: Contribution of neutrino and antineutrino modes compared to the total result assuming CPT 
invariance (solid curves) for is e disappearance (left panel) and disappearance (right panel) at the 90% 
CL, 2 d.o.f. . In both cases, 10 19 useful muon decays per polarity, = 2 GeV, and optimization point A 
( 20 m + 500 m) have been used. 



can easily see that the disappearance has a slightly better absolute performance, which 
mainly comes from the higher beam energy. In order to compare the absolute performance 
to existing experiments, we show the combined SciBooNE and MiniBooNE disappearance 



result from Ref. 25 , which was obtained in a similar spirit, as solid thin curve. The VLENF 
could improve this by about an other of magnitude. 

So far we have used the neutrino and antineutrino channels simultaneously, assuming CPT 
invariance. However, since the earlier mentioned anomalies appear for specific polarities, 
the separate sensitivities for neutrinos and antineutrino disappearance may be important 
for some models. We therefore show in Fig. [6] the contributions of the neutrino and an- 
tineutrino modes, together with the result assuming CPT invariance. As it can be clearly 
seen, similar limits can be obtained for the separated neutrino and antineutrino channels. 
The antineutrino sensitivities are somewhat weaker than the neutrino sensitivities due to 
the smaller cross sections. Nevertheless, CPT invariance (or other consistency) tests can 



be easily performed, see Ref. 13 for details. This is a significant advantage over many 
other experiments, which typically have the neutrino channels (beam experiments) or the 
antineutrino channels (reactor experiments) dominate. 

The discussed alternative setups only represent a limited selection of the ideas for sterile 



neutrino searches, see Appendix A of Ref. 11 for a more complete list. However, none of 



the proposed options seems to be able to compete with the proposed disappearance searches 
at the VLENF. 
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7 Summary and conclusions 



We have studied short-baseline electron and muon neutrino disappearance at a VLENF with 
Efj, = 2 GeV and = 4 GeV. Compared to the appearance channels, where backgrounds 
limit the sensitivities, the disappearance channels suffer from the knowledge on cross sections 
x efficiencies. We have therefore chose a setup similar to reactor experiments, such as Double 
Chooz and Daya Bay, using a near and far detector, in which the unknowns can be measured 
in a self-consistent way. For the systematics, we have adopted the most conservative case 
of a shape error fully uncorrelated among the bins, but fully correlated between the near 
and far detectors. In addition, we have included the extension of the decay straight and 
the beam divergence from the muon decay kinematics, which affect the sensitivity to very 
large Am 2 > 30 eV 2 in the near detector. Note that any additional machine-dependent 
divergence will add to the muon decay systematics and lead to some additional averaging. 
However, the muon decay kinematics cannot be eliminated, even if the other systematics 
might be improved. 

We have performed a two-baseline optimization of the setup, where we have identified 
optimal points depending on the value of Am 2 . From the different options, we have chosen 
a setup with the near detector as close as possible to the source and the far detector at 
a distance between about 500 and 800 meters, which is consistent with the optimization 



for appearance 12 and a good compromise between the small and large Am sensitivities 



As far as the minimal luminosity and muon energy are concerned, we have found that at 
least 10 19 useful muon decays per polarity are needed for E^ = 2 GeV, or, alternatively, 
a higher muon energy, in order to outperform practically any other proposed alternative 
setup and test the relevant parameter space. That is different from the appearance channel 
optimization, where lower luminosities may be sufficient and the muon energy hardly matters 



as long as E^ > 2 GeV 12 . Note that the VLENF setup can measure both electron 
and muon neutrino disappearance, for both neutrinos and antineutrinos. We have also 
demonstrated that the proposed setup is practically insensitive to the external knowledge 
on cross sections and efficiencies for Am 2 < 30 eV 2 , whereas the sensitivity for larger Am 2 
depends on the systematics assumptions since oscillations average out in both near and far 
detectors. 

In conclusion, it is well known from 3 + 1 models of sterile neutrinos that disappearance 
channels provide complementary and important information to constrain the models for 
sterile neutrinos. In the 3 + 1 case, the tension between the appearance and disappearance 
of various experiments, and between neutrino and antineutrino appearance in MiniBooNE 
have basically ruled out this model. Apart from the direct tests of disappearance anomalies, 
for more complicated 3 + iV scenarios, better disappearance information from both electron 
and muon neutrino disappearance will be needed. These channels may be uncorrelated in 
the underlying physics model, as it is evident already in the 3 + 1 case. The VLENF can 
provide this information if it is designed similar to the reactor experiments, with near and 
far detectors as similar as possible. 
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